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Abstract 


The  results  of  this  paper  concern  the  nonlinear  parabolic  P.D.E. 


-  5x  <VU)  -  £ 


u(0,  t)  =  u(l,  t)  =  0 
u(x,  0)  =  uq(x) 


where  u  represents  the  mean  velocity  field  of  a  turbulent  flow  in  a  plane  duct, 
and  vT(u)  =  G(x)  t)  |  is  a  model  of  turbulent  viscosity  based  on  the 

length  scale  of  this  flow. 

First  we  prove  the  existence  of  a  solution  as  well  as  same  results  on  re¬ 
gularity.  Then  we  demonstrate  that  u  is  positive  inside  the  domain,  hence  the 
same  result  holds  for  3u/3x(0,  t) .  The  latter  result  is  crucial  for  the  unique¬ 
ness  proof  which  is  developped  in  the  last  section. 
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A 


SIGNIFICANCE  AND  EXPLANATION 


_ \ 

The  turbulent  viscosity  assumption  is  the  simplest  turbulent  closure 
of  the  Navier-Stokes  equations.  Of  the  possible  nodels,  same  are  of  physical 
interest.  In  the  case  of  the  plane  duct,  a  model  based  on  the  length  scale, 
Uf,  of  the  flow  seems  appropriate.  Moreover,  the  introduction  of  such  a  model 
in  the  Navier-Stokes  equations  leads  to  a  novel  nonlinear  parabolic  equation, 
whose  mathematical  study  is  the  subject  of  this  paper. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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MATHEMATICAL  AND  NUMERICAL  ANALYSIS 
OF  A  TURBULENT  VISCOSITY  MODEL 

C.  M.  Brauner  and  C.  Laine 


i 
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0  -  INTRODUCTION 

This  paper  is  devoted  to  sane  ratharetical  and  numerical  results  on  a 
particular  modal  of  turbulence  in  a  plane  duct.  In  this  physical  case,  the 
statistical  average  of  the  basic  Navier-Stokes  equations  leads  to  a  one  di¬ 
mensional  mean  velocity  field  equation.  In  order  to  close  this  equation  we 
adopt  a  turbulent  viscosity  model  based  on  the  characteristic  velocity  scale 
Uf  =  (v0l3u/3xiwai1)  1//2  (see  the  Appendix  Al) .  This  efficient  model  induces 
a  nonlinear  parabolic  problsn 


(0.1) 


32u 


-  v  — 
3t  o  3x 


3x 


vT<u) 


G(x) 

u(0,t)  =  u(l,t) 


(VT(U)  S 
1/2 


0,  u(x,0)  =  uq(x) 


whose  study  is  the  purpose  of  this  article.  In  (0.1)  >  0  is  the  kinematic 

viscosity  of  the  fluid,  while  the  positive  function  G  adjusts  the  turbulence 
level  (see  details  in  the  Appendix) .  Several  assumptions  on  G  will  be  made 
below.  The  function  f  represents  the  pressure  gradient. 

In  section  1,  we  prove  the  existence  of  a  solution  via  a  Faedo-Galerkin 
method  and  we  construct  a  sequence  of  approximate  solutions  in  a  finite 
dimensional  space .  To  pass  to  the  limit  as  m  ■*  +  00 ,  we  make  use  of  a  compact¬ 
ness  argument.  In  section  2,  we  prove  further  results  of  regularity  by  the 
same  methods,  namely  in  the  space  L2  (0,  T  ;  H3  (f!))  n  L"  (0,  T  ;  H2  (fl) )  .  A  di¬ 
rect  consequence  is  that  u  is  continuous,  and  so  is  (0,  .) . 

The  main  point  of  the  paper  (Section  3)  is  to  prove  that  u  is  positive  in 
the  case  of  positive  data.  This  result  will  be  the  consequence  of  a  Hamack 
inequality  for  parabolic  equations  due  to  J.  Moser  [6],  [7] .  Then  we  adapt  a 
proof  by  Protter  and  Weinberger  [8]  to  demonstrate  that  (0,  t)  >  0.  The 

dX 

latter  result  is  crucial  for  the  proof  of  the  uniqueness,  which  is  developped 
in  the  last  section. 

Finally,  we  briefly  present  some  numerical  resuilts  in  Appendix  A2. 
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He  are  looking  for  a  solution  to  the  problem  : 


with  the  boundary  and  initial  conditions  : 


(1.2)  u(0,  t)  -  u(l,  t)  -  0 

(1.3)  u(x,  0)  -  uQ(x) 

In  (1.1),  vT(u)  is  the  “turbulent  viscosity"  : 

(1.4)  vT(u)  ■  G(x)  yigioTtTl 

let  n  ■  ]  0,  1  [  and  T  a  positive  nutber,  we  denote  by  Q  »  51  x  ]  0,  t[, 
and  we  assune  that  : 

(  f  6  L1  (Q) ,  u  €  H*  (£1) 

U-5J  {  ,  .  °  °  (+, 

(Gee1  (ff),  G(x)  >  Vj  >  o  1  ’ 

Let  us  introduce  the  following  notation  :  (  ,  )  is  the  inner  product  In 
L2  (£1)  and  1 1  .  |  |  the  associated  normj  |  .  1 1^  is  the  norm  in  L° (£1) .  Any  other 
norm  will  be  indicated  explicitly. 

We  introduce  the  eigenfunctions  of  the  operator  -  d2/dx2 ,  namely 
Wj(x)  ■  sin  jirx,  j  ■  1,  ...,  m  ...  .  Thai  -  w^"  *  X^wy  '*ere  X^  »  (jir)2. 

We  define  u^tt)  -  (t)w^  as  the  solution,  on  some  interval  [0,  tj, 

of  the  differential  system  : 

K(t)'  -  Vo  4r  <*>'  wj>  *  l&(0't)|1/2(|^G(x)  ^(t)V  Wj)  - 

(1.6) 

-  (f(t),  1(  j(m 


(+)  Practically  »  G(0)  «  v  ,  i.e.  near  the  wall  only  the  kinematic  visco¬ 
sity  occurs. 
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(1.7)  1^(0)  =  u0  =  Pr°jection  of  UQ  on  vm 

m 


<+) 


1.1  -  A  priori  estimates 


Lama  1.1  :  The  sequence  u^  is  bounded  independently  of  m  in  the  space 
L2  (0,  T  ;  tf  (ft))n  h  ( 0 ,  T  ;  L2  (SI)),  toreover  |j^(0,  .)(1/4  jp2  is  bounded  in 


the  space  L  (Q) . 


Ptoofj  Multiply  Eq.  (1.6)  by  g.  (t)  and  sum  from  j  =  1  to  m.  Then  ; 

■Tn 

.  .  3u  3u  ,  e  3u 

lit  *vo  Hsr1011’  *  isr10'  «i  Jti°<5rlt>,!  *•< 


where  c^  is  any  positive  constant.  According  to  Poincare's  inequality,  there 
exists  c_  >  0  such  that  v  1 1  g^lt)  |  | 1  5-  -p  1 1  u _(t)  |  | 2 


~2  1,1  H1  (fi) 

o 


Then,  choosing  =  c2^o'  an^  integrating  frcm  0  to  t  : 


ll^tJll2  t^pll  u  (t )  1 1 2  dr  + 

T"  c2  Jo  m  H1  (fi) 

o 

ft  3u_  ,  n  r  3u 

(1.9)  +2  |  (0,  t >  1  1/2  I  G  (r^(T))2  dxdr  ^  llujl2  + 

J  o  J  ft  in 

+  ^ffc  1 1  f  (T )  |  M  dr 


o  ‘o 

The  R.H.S.  of  (1.9)  is  bounded  by 
c„  iT 


ll«oll<+v“ 

n 


1  f  (t>  1 1 2  dt  =  c. 


o  ‘ o 


As  1 1 u.m (t)  1 1 2  ^  c^  V  t  e  [0,  T  i,  we  infer  that  t^  =  T  and  we  have  the  estimate 
in  L“  (0,  T  ;  L2  (fi)) .  Talking  t  =  T  in  (1.9)  gives  the  estimate  in  L2  (0,  T  ;  H^(fi)) 
Finally  as  G  >  v1  >  0  : 


(+) 


Vm  is  the  space  generated  by  Wj^ 


,  •  •  •  r 


(1.10)  JQ  1 5?  (0,  t)  1 1/2  (g^  (t))2  dx  dt  «  ^ 
and  the  proof  of  the  lama  Is  conplete. 


Now  we  establish  a  sharper  estimate  : 


Lama  1.2  :  Hie  sequence  u^  is  bounded  independently  of  m  in  the  space 
L2(0,  T  t  H 2  (£1) )  r>  L"(0,  T  ;  Moreover  <0,  .)1/4  |—  is  bounded 

in  the  space  L2  (Q) . 


Pk  oo  A  :  In  (1.6)  we  replace  w..  by  -  1/Xj  w.."  ;  we  multiply  (1.6)  by  g.  (t) 
and  again  sum  fran  j  =  1  to  m.  Clearly  : 


(1.11) 


A  3u„  3  Uu,  .  3u  ,,,  ,  3  u 

ft  ll|?  <t)||2  +  2vq  1!^  (t)||  +  2  (0,  t)|1/2  J^G(^(t))^  dx 

3u  1  n  t  Af  3u  a2u  *  2 

+  Iif(0'  «l1/2-j niar(t)  jsr™  dx  =  2(f(t),  y»(t)) 


Integrating  fran  0  to  T,  it  turns  out  that  : 

3u  32u_  3u 

II|?W|I^2v0||^(21(q)+2v1  ||  |^<0,.> 

3  2u 

(1.12)  <  2  l|f||T2,ft.-  11-3=1 


1/4  a*um,,t 


33?^"L2(Q)  4 


l2(b,|IS?=iii,(q)  +  HeBT211  + 


3u 


32u 


ln,(0,t)|1/2  |g|  |J?(t>|  |^(t)|  dx  dt 


Wfe  can  find  an  upper  bound  for  the  R.H.S.  in  the  following  way  : 

3 vi  .  y_  „  3u  32u 

1^(0,  t)  | 1/2  |£|  \^t)\ 

JQ 


T  9 vi  .  ^  3u  3zu 

x  -  2  ]  lj?(o,  t) | 1/2  |g|  |^(t)|  dxdt  = 

(1.13)  (  -  2  J  |g  [lj?(0,  t)|1/4  |J?(t)||  t)|1/4  |^(t>"j] 


dx  dt 


<  2  HflL  M  l5?(0-)l1/45?l 

and  by  Young's  inequality. 


3u 


2  II  I5~(o,.)|^3-1hii 

L2  (Q)  3x  3x  L2  (Q) 


,  /i.  3*u 
1 1/4  mi 


-4- 


+  -  1 1—1 
v1l 'dx' 


m,n  v  |  1/h  _ mi 

— (°,-)|  'l2  (q) 


c»U  i  /A  r)U  2 


Now  taking  (1.14)  into  account  and  again  applying  Young's  inequality  to  the 
quantity  1 1 f  |  lL2  (Q)  .  1  lip21!  1 L*  (Q) »  U  -12)  beccmes  : 


Ill55f(T)||2  +  vo|l^HL2(Q)  + 

+  V1  II  I  33r(0,t)  I  ^  (ax5-  ^  1  ^L2  (Q)  *  vo  I  ^L2  (Q)  + 

As  u_  -►  u  in  H' (fi) ,  by  (1 .10) , 
cm  o  o 

3\ 

(1.16)  .-7s1  is  bounded  In  L2  (Q) . 


Now  the  estimate  in  L°(0,  T  ;  H1  (ft))  is  obtained  by  integrating,  (1.11) 

O  gj.  1/4  3  Um  2 

fran  0  to  t,  0  <  t  <  T.  Finally,  let  us  note  that  11  | ^(0, . )  |  '  ■  |  2 , 


is  bounded  by  (1.15)  ;  so  Laima  1.2  is  proved. 


3x^'  'L2  (Q) 


Lama.  1 .3  :  The  sequence  u  '  is  bounded  in  L2  (0,  T  ;  H  1  (fi) ) . 

-  m 

Paopj  i  Introducing  the  projector  Pm  fran  L2  (fi)  into  V  ,  we  may  write  : 


ir10'  l)|1/2k<GS!)]*pn£ 


It  is  well  known  that  Pm  enjoys  the  following  properties  :  1 1  Pm|  Lg(Hi  .  Hij 
and'  as  Pm  is  ^f^oint,  l|Pmll^(H-l  .  H"l)  o  o 


As  a  result  of  Lama  1.2,  u^  is  bounded  in  L2  (0,  T  ;  Hl(fi))  and  by  Sobolev 
embedding,  in  L2  (0,  T  ;  C1  (R) ) .  Therefore,  ISiy^xfO, .)  | 1//2  is  bounded  in  L*(0,T) 
then,  as  (G  3u  /3x)  is  bounded  in  L°(0,  T  ;  H_1(fi)) : 


(1.18) 


,av 


1^(0,.) 


|l/2  L.  (a  ^ 

’  3x  3x  ' 


Is  bounded  in  L4(0,T  ;  H-1(ft)), 


hence,  Lenina  1.3  is  proved. 

■ 

1.2  -  Passage  to  the  limit 

Rom  the  above  lemnas,  it  follows  that  we  may  extract  a  subsequence 
such  that 

(1.19)  +  u  in  L2  (0,  T  ;  H2(0)  fl  H^(n))  weakly. 

By  a  oonpactness  theorem  (see  e.g.  Lions  [4],  p.  57-58  with 
BQ  -  H2  (£)) ,  Bx  -  H“l(£)),  B  -  H2"6^),  e  >  0), 

we  can  extract  another  subsequence,  also  denoted  by  u^,  such  that 

(1.20)  -*•  u  in  L2  (0,  T  }  H2-6^))  strongly. 

But  we  may  choose  e  >  0  small  enough  that  the  Sobolev  embedding  yields 
H2-6^)  C  Cl  (S),  (in  facto  <  e  <  |).  Then 

(1.21)  -»  u  in  L2  (0,  T  ;  c1^))  strongly 

and 

(1.22)  1^/(0,. )|1/2  -*  ||H(0,.)  1 1/2  in  L*  (0,T)  strongly. 

Let  j  be  fixed  in  (1.6)  with  j  <  y.  As  y  we  get 

(u'  (t) ,  Wj)  -  vo  (^(t),  w.)  -  ||H(0,  t)|I/2  (|j  (G  §jj(t»,  Wj)  -  (f(t),  w.) 

for  j  *  1,  2  ...»  hence  u  satisfies  the  equation  (1.1). 

So,  we  have  the  following  existence  theorem  > 

Tke.o >iem  1.1  j  under  Assuiptions  (1.4)  and  (1.5),  Problen  (1.1)  ...  (1.3) 
possesses  a  solution  in  the  space  L2(0,  T  ;  H2(fl))  D  l"(0,  T  ;  H^(0)).  Moreover 

||H{0, .)  | 1/4  and  |||(0, .)  | 1/4  |jl  belong  to  L2  (Q) . 


-6- 


In  order  to  prove  the  uniqueness,  further  regularity  results  nust  be  es¬ 
tablished. 


2  -  RESULTS  CM  REGUIARTIY 

2.1  -  Regularity  in  Sobolev  spaces 

Th.zon.em  2.1  :  Let  us  assume  that  f ,  G,  and  uQ  are  given  functions  with 

(2.1)  f  e  L2 (0,  T  ;  H^(ft) ) 

(2.2)  G  6  C2(S3)  ;  G(x)  »  v]L  >0,  g'(0)  =  G’(1)  =  0 

(2.3)  u  eH2(«)nH'(fi) 

o  o 

Then  any  solution  u  to  problem  (1.1)  ...  (1.3)  obtained  via  the  preceding 
section  satisfies 

(2.4)  u  6  L2(0,  T  ;  B3(fi))  D  L°(0,  T  ;  H2 (fl) ) 


-7- 


The  main  term  is 


-  ( ^  (e  ^t»,  ^(t»  =  -  ( ^t)) 


% 


Du 


3  u 


dx  3x 


3x4 


32u  3  u 

(  G  3P"lt)'  3F“(t))' 


we  compute 


1  dx  Dx  lt) ' 


3  u_ 
33? 


J(t)) 


»  Du  ,  3s u 


dx  Dx 
as  G1  (0) 


3x3 
G*  (1) 


.2„  3u  3’u  ^  3’u  33u 

<  3#K>’  a?®'")  *  <  Si?®'11'  Jp5'*'’- 


and. 


32u  3“u  a  D2  u  33u 

-  <  G  3^(t)'  3^(t,)  *  (  h  (G  3F2-(t,)'  3^(t)) 

„  32u  33u  33u  33u 

=  <f  5^(t) '  3^(t)>  +  (G  3^(t> '  S*T(t)> 


(2.7) 


1_  d_ 
2  dt 


m(t)||  + 


Consequently  (2.5)  may  be  written 

32u  2  3Ju_ 

I^H  +vo 
3u  1/2  D’u  3Ju 

+  I  10»  fcM  (  G  1 3x1“^))  : 

s  -  2  |^p(0,  t)|1//2  <dG  '^n' 


_  |!^n(0f  t)  | 1/^2  '  d2G'*Um,*.»  Un> 


<*?*?<«•  st®"=»  -<Slu>5^!|t» 


By  integration  from  0  to  T,  it  follows  that 


(2.8) 


D2u  * 

+2vo 


33u 
I.  IS 
Dx5^ 


3u 


+  2 


|  |^(0,  t)|1/2  G  (  ^(t))'  dx  dt  = 


'L2  (Q) 

33u 

m. 


¥ 

D2u_  2  I  3u  3su  j3u 

=  11^0)11  -  4  (Q  !af(0'  fc)l  S  ^(t) 

9u 


-  2  <**  -  2  j  |^(0,  t)\V-^  g^(t)  5^(t)  dx  dt 


-8- 


Assumption  (2.2)  about  G  then  enables  us  to  get  : 

/ 1  N 2  +  2^o  1 'l2  r0)2  1 1  '&(o'  ■ •>  i 1/4  S^i  *l2  (q)  « 

\  .<11^(0,11%  4  HfiL  II  !^(0,.)|1/4^I!l2(q)  . 


1 3x  <0'*)l  ^  3x3  I  I L2  (Q)  +2  I  I  3x1  I L2  (Q)  ‘  I  1 3x3  ^  '  L2  (Q)  + 


1  dx  00  3x 


3u  ,  ..  3u„  3u  .  33u 

Sx-*0'^!  3x~ I  ‘L2  (Q)  II  I  3jc^(0,  I  3x^1  I L2  (Q 


By  applying  Young's  inequality  to  the  R.H.S.,  we  have,  with  c^being  any  positive 
constant , 


3u  ,  3  u  3u 

"liui  i5J!<o.->i‘/)3?!!i:L.ra  n  W 


<  i^iSnii  ig^r 


3u  ,  ..  33u 

^T(0'-)l  3x  “I  I L2  (Q)  * 


3x3  1  1  L2  (Q) 


+  2ciiiiiin  i^^i^^'^q, 


On  the  same  way,  for  any  > 


dx3  *  oo  3x 


<  iiiUnii  iw5<».-n‘/,d!ii:!,„,  * 


3u  .  3u  3u  ,  ..  33u 

ll}f1<0'‘)l  3x~l  II2  (Q)  'I  lsjf1(0',)l  3x^1  I L2  (Q)  ^ 


C2!  'dx2  lloe  3x 


i  ,  ,A  ou  2 

5i(o  \  i J/4 _ nil  i  , 

.  l0'*>  I  -dx  I  lL2  (Q) 


“2'  'dx3  L  1 1  3x 


and  for  any  c3  >  0 


m,n  %  1 1/4  Um,  | 2 

a^ML2(Q) 


2  I  1 3x I  I L2  (Q)  I  1 3x3  I  II2  (Q)  ^  c3  I  Isxl 'l2  (Q)  +  C3  I  Isx^l  II2  (Q) 
Then  it  is  convenient  to  choose  Cj ,  and  c3  such  that 


c3  =  vo  '  2ci  i^ll  +  c: 


2  "d^L  =VI 


Fran  the  theorem  of  intermediate  derivatives  (LI0NS-MW3EUES  [5])  we  infer 
that  u  belongs  to  c°([0,  T]  ;  H2(ft))  and  by  Sobolev  embedding 

(2.12)  u  6  C°([0,  T]  ;  C*(Q)) 
hence  the  theorem. 


3  -  MMCMM  PRINCIPLE 

In  order  to  prove  uniqueness,  it  is  crucial  to  establish  the  positivity  of 
3u/3x(0, .) .  This  result  will  be  the  consequence  of  the  positivity  of  u  in  Q. 

This  latter  result  oanes  fran  a  Hamack  inequality  for  parabolic  equation  due  to 
J.  MOSER. 


3.1  -  A  Hamack  inequality 

Let  us  recall  the  following  result  due  to  J.  MOSER  [6], [7]. 

Let  n  denote  an  open  bounded  set  in  IRn  and  Q  =  !lx  ]0,  T[.  We  consider 
the  equation  in  Q 

(3.1)  -  l  (aij  (x'  t}  3^  ^  =  0 

J-r  J  J  J- 

where  a^  are  bounded  functions  such  that  : 

°  <  X  <  I  aij  Ci  »  V  e  =  (C1#  ...  5n)  e ,  lei*  -  1 

j 

(3.2) 

aij  =  aji 

Let  u  be  a  vreak  solution  of  (3.1)  such  that  3u/3t,  Su/Sx^  ,  i=l,...,n  , 
belong  to  L2 (Q) .  We  suppose  that  u  is  non-negative  in  Q.  We  consider  a  caipact 
and  connected  subdomain  of  Q,  and  two  subintervals 

i”  =  (t,  tj  <  t  <  tj)  ,  I+  =  (t,  t3  <  t  <  t4)  where  we  assume  that 
0  <  tj  <  t2  <  t3  <  t^  <  T.  Define 

Q"  =  I-  x  A  ,  Q+  =  I+  x  A 
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lama  3.1  :  Any  non-negative  solution  u  of  (3.1)  satisfies 


(3.3)  ess  sup  u  <  c  ^  +  *  *  ess  inf  u 

Q"  Q+ 

where  c  >  1  is  a  constant  which  depends  only  on  Q,  Q+,  Q  . 


3.2  -  Positivity  of  u 

Here  ft  is  again  1 0,1[ .  We  shall  apply  the  above  result  of  J.  MOSER  to  an 

auxiliary  problon.  First  we  prove  that  u  is  non -negative  v-ta  the  weak  maximum 
principle. 

lama  3.?  :  The  hypotheses  are  (1.4),  (1.5).  Besides,  suppose  ^non-nega¬ 
tive  in  ft  and  f  non-negative  a.e.  in  Q,  then  any  solution  u  of  Problem  (1.1)  ... 
(1.3)  is  rcn -negative  a.e.  in  Q. 


Paaoj  :  As  in  LIONS  [4]  p.  290,  we  multiply  (1.1)  by  u~(t).  Then  after 
Integration  over  ft,  it  follows  that  : 

-It  I  l«~<t)l  |2  -  2*0  llf£  (t)H2  +  2  Jn  l|^<0,t)l 1/2  G(|a  tt))*  dx 

-*■  2  f  f  (t)  u"  (t)  dx 

Jn 

which  yields  ^  |  |u_(t)|  |  <  0.  to  u’  «  0,  1 1  u—  (t)  I  I  ■  0,  V  t. 

■ 

For  u  any  solution  of  Problem  (1.1)  ...  (1.3),  we  introduce  now  the  auxi¬ 
liary  linear  problem 


l  u(x,  0)  -  uc(x) 


Luma  3,3  :  The  hypotheses  are  (2.1),  (2.2),  (2.3),  u  given  by  Theorem 
2.1.  Also,  suppose  f  non-negative  a.e.  in  Q.  Then  u  e  C® (Q)  and  u(x,  t)  >  u(x,  t) 
v  (x,  t)  e  Q. 

Ptooi  s  Problem  (3.4)  admits  a  unique  solution  u  e  L2  (0,  T  ;  HJ  (ft)  O  H1  (ft)) 
•v,  -  o 

!£e  L*(Q),  then  u  e  C®  (Q) . 

-12- 


Set  w  =  u  -  u  and  w  satisfies  : 


(3.5) 


a?  ‘  vo  &  ~  (i(0- 1))1/2  h  (G(X)  "  f  *  0 

w(0,  t)  =  w(l,  t)  =  0 
w(x,  0)  =  0 


Clearly  w  >  0  by  the  v«ak  maxinum  principle. 


Thexixem  3. 1  :  The  hypotheses  are  (2.1),  (2.2),  (2.3),  u  given  by  Theo¬ 
rem  2.1.  Besides,  suppose  f  non  negative  a.e.  in  Q,  and  uq  positive  in  ft.  Then 
u  positive  in  Q. 

:  By  Lama  3.3  it  is  sufficient  to  prove  that  u  >  0  in  Q.  We  apply 
Lanma  3.1  to  u  with  ft  =  ]n,  i[,  a^  (x,  t)  =  a(x,  t)  =  G(x)  (i^(0,  t))1/2  +  vq, 

A  =  v  and  A  *  v  +  {  sup  G(x)  .  sup  (|^(0,  t))1/2  l. 

°  x  e  ft  t  6  [0,  T]  3x  ’ 

We  prove  u  >  0  in  Q,  exactly  as  we  did  in  Lenina  3 .2 . 

New  let  xq  e  ft  be  fixed.  Since  u  e  C°  (Q)  and  uq  >  0  in  ft,  there  exists  a 
closed  interval  A  c  ft  containing  xq,  and  a  time  T,  such  that 

(3.6)  u(x,  t)  >  0  V  x  e  A,  O^t^T 

Next,  let  t#>  x  be  such  that  u(xQ,  t*)  =  0.  With  the  notations  of  Laima 
3.1,  we  choose  ^  and  such  that  0  <  tx  <  t2  <  t.  Then  we  consider 
t3  *  t#-  e,  t^  *  t#+  e  (e  >  0  anall  enough) . 

So,  infQ+  u  =  0.  Therefore  as  a  result  of  (3.3)  supQ_  u  =  0  which  is  ireieed 
inconsistent  with  (3.6) ,  and  there  does  not  exist  such  a  time  t#. 


Rematfe  3 . 1  :  The  above  demonstration  may  be  extended  to  the  caseu  >  0  only 
on  a  subdomain  of  ft. 

3.3  -  Bositiyity  of  3u/3x(0..) 

Fran  the  positivity  of  u,  we  will  deduce  the  positivity  of  3u/3x(0, .) . 


We  now  Introduce  the  auxiliary  function  v(x,  t): 
v(x,  t)  =  -  exp  (-8  [  (x  -  Xj)2  +  (t  -  t^) 2 1 )  +  exp  (-  8  R2)  and  note  that 
L  v  =  28  exp  (-  8[(x  -  x^2  +  (t  -  t^2])  [28a  (x  -  Xj)2  -  a  -  b(x  -  x^  +  t -  t^ 

Thus,  for  sufficiently  large  8  we  have  L  v  >  0  in  D.  Set  w  =  u  +  e  v  for 
£  >  0.  Thus  Lw=f  +  eT,v>0  a.e.  in  D.  Because  of  fact  (iii) ,  we  can  choose 
e  so  stall  that  w  >  0  on  C'.  Since  v  =  0  on  3K,  we  have,  because  of  (i)  w  >0  on 
c”»  except  at  P  and  w(P)  -  £> 

It  is  easy  to  shew  by  the  weak  naxiitun  principle  that  w>  0  In  D.  We  deduce 
that  |^(P)  5.  0. 

On  the  other  hand,  we  see  that  (3v/3x) (P)  <  0  since  : 

(x,  t)  =  28  (x  -  xx)  exp  (-  e[(x  -  x^2  +  (t  -  tx)  ), 

§£  (P)  =  -  28  x:  exp  (-  8  R2)  <  0,  and  finally,  |£(P)  =  |£(P)  -  e  |X(P)  >  0 

Co>io£Jta/iy  3.1  :  the  hypotheses  are  those  of  Theorem  3.2.  Moreover, 
suppose  is  positive.  Then  there  exists  a  constant  a  >  0  such  that 

(3.7)  |^(0,  t)>a>0  V  t  6  [0,  T] 

Ptoo((  :  The  corollary  is  an  easy  consequence  of  Theorem  3.2  and  of  the 
continuity  of  (0,  .) . 


4  -  UNIQUENESS 

The  above  results  will  enable  us  to  prove  the  uniqueness  of  a  solution 
to  Roblan  (1.1)  ...  (1.3),  namely  : 

Theorem  4.1  :  The  hypotheses  cure  (2.1),  (2.2),  (2.3) .  Besides  we  assure 
that  f  is  non-negative  a.e.  in  Q,  uQ  is  positive  in  fi  and  g£P(0)  is  positive  , 
Then  the  solution  of  (1.1)  ...  (1.3)  is  unique. 

PKQoj  :  The  proof  of  the  uniqueness  is  a  non  standard  application  of 
Gronwall's  lama. 
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Suppose  there  exist  tvc  solutions  and  U2 •  Summarizing  the  results  of 
section  3,  there  exist  2  positive  constants  a  arxi  K  such  that  : 

(4.1)  0  <  a  $  (|j*i<0,t))1/2  $  K,  1=1,2  V  t  8  [0,T] 

Now  set  u  =  Uj  -  Uj,  and  u  satisfies 

(4*2)  + 

+  {^  <0'  «)1/2  "  ^  (0,  t))1/2}^  (G  |f) 

(4.3)  u(0,  t)  ■  u(l,  t)  =  0 


(4.4) 


u(x,  0)  =  0 


HUtiplying  (4.2)  by  -  ~  yields  s 


(4.5) 


is  llgwll*  *  »o  <g|£,  |i)  - 

-  (^(0,  t))W  (*  |M,  . 

-  (gMo.  t))1-*  -  (|Hi(0,  „)«)  (f  |i“) 

-  <C|gM0.  t:)1^  -  <|S2(a,  ti)1'2)  (g§2,  g?) 


Note  that  : 


(4.6) 


l<&».  t))1/2  -  <ff(0,  l?x{0'  t}  -5^(0>  fc)l 

(|^(0,  t))1^  (|f(0,  t))1/2 


.  j£gl  t}  ||f),  t)| 

a  'a 

Since 

«•’>  ifi«.wu  V7ii!!<»ii1/2ii&<t>ii1'2. 

it  follcws  that 
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|(|HI(°,  t))1/2  -  (|H2(0,  t))1/2|«  £  ||g(t)||1/2  ||0(t) 


and  by  Insertion  into  (4.5)  : 


1  d  1 1 3u 

2  dt  ' '5x 


(fc)  1 1 2  +vQ  llf^<t)||2  +a  ||0(t)||2  < 


*  £  I  @1-  I  ^(t)  1 1  !  1 11/2  1  & lV2 

,  IlGlI.  ||gftt>||  ||g(t)||l/2  ||0(t)||3/2 

To  estimate  the  R.H.S.  we  apply  the  Young  inequality  three  times  :  C1#  C2,  C3 
being  as  usual  any  positive  constants, 


. . * 

♦¥"§"-  notin' 


•  <  ret  *  Sjllljfltlll  *  rf-  ll$g)(t)ll 


4  C2C3  1  1 3x 


Therefore  (4.8)  becomes,  since  ^2  e  L° (0,  T  ;  L2(fi))  from  Theorem  2.1, 

I  i-ft  -  n|^<«n'[v0-«,  v-i  ugii_ - 

l  M®"- 1  *' 1  % *  ^ 

I  <  US*011  [^7  1  '^1  I-  +  {!  ,G|  1  'l^1  ^-(o.T.Ll 

\  *  Hill-  IISPIU 


♦  n§iuil?i 


L  (0,T;L2  (fi) )' 


Choosing  C3 ,  C 2 ,  C3  in  order  to  eliminate  the  coefficient  of  1 1  ^  (t)  1 1  in  the 
L.H.S.,  it  follows 


Srllgtt>ir<<*  I ISWII 


dt  11 3x 


Plan  Grcnwall’s  lama,  |  ||£(t)|  j  =0  ,  Vt.  Hence  |jjl  =  |Jj2  in  Q,  ani 
the  uniqueness  is  proved. 


APPENDIX 


A1  -  THE  PHYSICAL  PR0BUEM 

Hie  flow  under  consideration  is  incompressible  ;  the  equations  describing 
the  turbuloit  field  follow  from  tin  statistical  iver age  of  the  Navier-Stokes 
equations,  using  the  classical  decomposition  of  the  velocity  field  into  mean 
and  fluctuating  components  =*  0^  ...  Mb  have  then  : 

30,  _  3U,  *  320,  _ 

<,v*1,  at-  +  Uj  "  p  3xt  +  vo  3xi*'  “  3Xj  (uiuj) 


To  close  Equations  (A.l)  and  (A. 2) ,  we  use  the  turbulent  viscosity  asstnp- 
tions  given  by  : 

_  30,  30,  .  _ 

CA.3)  -  ulUj  .  vT  (^.  +  jjl)  +  j  q2  «tj 

where  u^Uj  is  the  Reynolds  stresses  tensor.  In  (A.l)  the  kinetic  energy  q2  =  u^ 
will  be  ahecjrbed  into  P  and  so  will  not  need  to  be  calculated  explicitly. 

In  particular,  we  consider  the  turbulent  fluid  in  notion  between  two  paral¬ 
lel  planes  at  a  distance  2  D.  Let  ufxj,  t)  be  the  mean  velocity  of  the  fluid  at 
a  point  6  [0,2Dl ,  at  a  time  t.  The  physical  simplifications  of  the  general  e- 
quations  lead  to  : 

..  3U  32U  1  3P  3  , - , 

<A,4)  3t  ’  vo  3jF  *  ”  p  3x:  “  3Xj  <U1U2) 


In  this  case,  the  turbulent  viscosity  \>T  is  chosen  as  a  function  of  the 
characteristic  scale  uf,  and  so  is  defined  by, 

«*•* 


where  G(x)_ls  a  polynomial  function  adjusting  the  turbulence  level.  Consequently, 

1  3P 

with  -  —  ^  »  f  being  given,  and  using  the  turbulent  viscosity  model  defined  in 
(A. 5) ,  we  hAve 

<*•« 

where  x  henceforth  only  denotes  the  spatial  variable  x^  »  which  is  here  a  dimen¬ 
sionless  variable. 


A2  -  9CME  NUMERICAL  RESULTS 

We  are  interested  in  confuting  the  established  solution  of  Problem  (1.1) ... 
(1.3)  : 


(A.l) 


~  v°  5yy  =  *G<y)  ! Cy  <°>  1 1/2  Uy)y  +  f  on  Si  =  ]o,  2D[ 
0(0)  =  U(2D)  =  0 


First  we  reduce  (A.l)  to  dimensionless  equations  by  introducing  the  cha¬ 
racteristic  length  and  velocity  of  the  flow,  defined  by  : 


(A. 2) 


L  =  2D 
o 


uf  myi 


dx 


Using  a  new  set  of  dimensionless  variables. 


(A.3)  u  -  0  /  uQ  ,  x  =  y/^  ,  Rf  =  (uf.2D)/vo, 

the  mean  velocity  field  beocnes, 

(  -  uxx  =  Rf1/2  <G<*>  I  V0)  1 1/2  Vx  +  2  Rf  011  ]0'  ^ 

l  u(0)  =  u(l)  =  0 

in  which  we  choose  G(x)  according  to  a  modified  Van  Driest  model  ([9]  p.  194) 

The  computation  has  been  realized  with  a  classical  fixed  point  method  ; 
also,  we  have  tested  an  optimal  control  method,  to  point  out  the  capability  of 
this  algorithm  to  solve  nonlinear  equations  with  a  very  large  nonlinear /linear 
ratio  (in  the  experimental  cases  Rf  is  larger  than  5000  1 ) . 
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We  have  used  a  non  regular  siixUvision  of  ]0,1[  In  order  to  have  a  better 
approximation  of  the  variables  in  the  wall  regions ,  where  the  gradients  are 
large.  Namely  we  take  ax2  -  subdivision.  The  Internal  approximation  la  reali¬ 
zed  with  the  Lagrange  finite  element  PI .  The  algorithm  efficiency  depends  on 
the  physical  meaning  of  initialization .  The  "parabolic"  profile  corresponding 
to  the  laminar  solution  leads  to  good  convergence  of  the  method  (note  that  the 
sequence  Is  alternating)  , 

Physically,  this  model  has  a  good  efficiency  since  it  does  not  demand  any 
parameter  adjustment.  Moreover,  it  Is  able  to  predict  the  mean  velocity  field 
Into  the  whole  duct.  Other  tested  models  only  work  In  a  subregion  of  the  duct, 
whose  border  belongs  to  the  so-called  "logarithmic  zone”.  In  this  case,  lnhomo- 
genous  boundary  conditions  should  be  prescribed.  Figure  2  shows  a  comparison 
between  the  experimental  results  of  OCHTB-BEUOT  [1]  and  our  numerical  results 
for  a  Reynolds  lumber  Rf  »  5160. 


W2EWJ32S55  -  This  paper  was  conpleted  when  the  first  author  was  visiting  the 
Mathematical  Research  Center  of  the  university  of  Wisconsin-Mad  Ison.  We  are  indtb- 
fcad  to  Professor  M.G.  Crandall  who  brought  J.  Moser's  work  to  our  attention. 


[1]  OCMIS-BELLOT,  G. ,  Ecoulanent  turbulent  entre  deux  parols  planes  pared. - 
Idles,  Pub.  Sci.  Tech,  du  Min.  de  l’Air,  419  (1965). 

[2]  LAINE,  C. ,  Etude  mathdmat ique  et  nunSrique  de  noddles  de  turbulence  en 
conduite  plane,  Tbdse,  University  de  Lyon  (1930). 

[3]  LAINE,  C. ,  Sur  les  aspects  mathdratlques  de  noddles  de  vlscositd  de  la 
turbulence  en  Seoul  anent  Stabli  de  conduite  plane.  C.R.  Acad.  Sc.  Paris, 
sdrie  B,  291  (1980) . 

[4]  LICKS,  J.L.,  Quelques  mdthodes  de  resolution  de  prcbldmes  aux  llmites 
non  lindaires,  Dunod  (1968) . 

[5]  LIONS,  J.L.  and  E.  MAGENES,  Prcbl&nes  aux  llmites  non  homogdnes  et  appli¬ 
cations,  Dunod  (1968) . 

[6]  MOSER,  J. ,  A  Hamack  inequality  for  parabolic  differential  equations, 
CCnm.  Pure  Appl.  Math.,  17,  (1964)  pp.  101-134,  and  correction  in  20 
(1967),  pp.  231-236. 

[7]  MOSER,  J.,  On  a  pointwise  estimate  for  parabolic  differential  equations. 
Conn.  Pure  Appl.  Math.,  24  (1971),  pp.  727-740. 

[8]  PROCTER,  M.H.  and  WEINBERHR  H.F.,  Maxinun  principles  in  differential 
equations,  Prentice-Hall,  (1967). 

[9]  TAHDR,  C.  ;  HUMS,  T.G.,  MDRGAN,  K.,  A  miner ical  analysis  of  turbulent 
flow  in  pipes,  Corputers  and  Fluids,  vol.  5,  pp.  191-204,  Pergamon  Press, 
(1977) . 


■22' 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  DM e  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

2366  f; l-fr/U 

S.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (md  Subtitle) 

MATHEMATICAL  AND  NUMERICAL  ANALYSIS  OF  A 

TURBULENT  VISCOSITY  MODEL 

S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 

6-  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORf.J 

C.  M.  Brauner  and  C.  Laine 

•.  CONTRACT  OR  GRANT  NUMSEROJ 

DAAG29-80-C-00  41 

S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 

610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 

to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

1  -  Applied  Analysis 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 

P.O.  Box  12211 

Research  Triangle  Park,  North  Carolina  27709 

12.  REPORT  date 

April  1982 

is.  number  of  pages 

22 

14.  MONITORING  AGENCY  NAME  A  ADDRESS (It  dtltoront  horn  Controlltnd  Ottlco) 

IS.  SECURITY  CLASS,  (ot  thle  report) 

UNCLASSIFIED 

1  So.  DECL  ASSI  FI  C  ATI  ON/  DOWN  GR  ADIn  G 
SCHEDULE 

1 16.  DISTRIBUTION  STATEMENT  (of  thle  Report)  1 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  tfio  ebetract  entere d  In  Block  20,  If  different  from  Report) 


IS.  SUPPLEMENT  ARY  NOTES 


It.  KEY  WORDS  (Continue  on  teweree  aide  If  neceeeary  and  Identify  by  block  ntenber) 

Turbulent  viscosity  model,  nonlinear  parabolic  P.D.E.,  maxinun  principle. 


20.  ABSTRACT  (Couth mm  on  tor  moo  tldo  It  n.  ceeemr  *nd  identity  by  block  number) 

The  results  of  this  paper  concern  the  nonlinear  parabolic  P.D.E. 


du 


o  U  3  /  ,  ,  3U\ 

3t  ~  Vo  to?  "  3x  ^T(u)  3x^ 


U(0,  t)  ■  u(l,  t)  *•  0 
u(x,  0)  »  u  (x) 


DO  , 


FOAM 

JAN  7S 


U71  EDITION  or  I  NOV  St  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whan  DMe  Entered) 


1 


vAere  u  represents  the  mean  velocity  field  of  a  turbulent  flow  in  a  Vna  due. 
a™*  v^(u)  ■  G(x)  y  lg“(0,  t)|  is  a  rnodel  of  turbulent  viscosity  basal  on  the 
length  scale  of  this  flow. 


^rst  we  prove  the  existence  of  a  solution  as  well  as  sene  results  on  re¬ 
gularity.  Then  we  demonstrate  that  u  is  positive  inside  the  danain,  hence  the 
W®  result  holds  for  3u/3x(0,  t) .  The  latter  result  is  crucial  for  the  unique¬ 
ness  proof  which  is  developped  in  the  last  section. 


